Effects of Local Fields on Spontaneous Emission in Dielectric Media 
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The local-field renormalization of the spontaneous emis- 
sion rate in a dielectric is explicitly obtained from a fully mi- 
croscopic quantum-electrodynamical, many-body derivation 
of Langevin-Bloch operator equations for two-level atoms 
embedded in an absorptive and dispersive, linear dielectric 
host. We find that the dielectric local-field enhancement of 
the spontaneous emission rate is smaller than indicated by 
previous studies. 



In the formative period of nonlinear optics, Bloember- 
gen taught us that the nonlinear optical effects of a dilute 
collection of atoms that are embedded in a dielectric host 
are enhanced by local-field effects jjj . Now, in the era of 
quantum optics, researchers are again looking at the in- 
teraction of dielectric materials, the radiation field, and 
resonant atoms. Central to these investigations are ef- 
forts to quantize the electromagnetic field in dielectrics. 
One widely used technique is to quantize the macroscopic 
Maxwell equations in which the classical constitutive re- 
lations have been assumed In the microscopic 
approach, a generalized Hopfield transformation, based 
on Fano diagonalization, is used to obtain the coupled 
polariton modes of the coupled field-oscillator system 
|^| flO|j . These quantization methods are well-established 
for dielectrics with negligible absorption and techniques 
to deal with the special requirements of quantizing the 
electromagnetic field in absorbing dielectrics are begin- 
ning to emerge |3]-@,[) 10 1. 

Since Purcell first predicted the alteration of the emis- 
sion rate of an excited atom due to an optical cavity [pd| , 
it has become well known that the observed spontaneous 
emission rate of an atom depends on its environment. 
When the quantized coupled field-dielectric theories are 
applied to the spontaneous emission of a two-level atom 
embedded in an absorptionless dielectric, the relation 
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is obtained P,^,[7|P JT(i|j . Here, n is the linear index 
of refraction and £ is the dielectric local-field enhance- 
ment factor, Tq is the vacuum spontaneous emission rate, 
and Tgg is the enhanced spontaneous emission rate in 
the dielectric. Both the Lorentz virtual-cavity model 
£ = (n 2 + 2)/3 and the Onsanger real-cavity model 
£ = 3n 2 /(2n 2 + 1) have been utilized in various studies 
of local-field effects on spontaneous emission. 

One of the key features of these approaches of apply- 



ing the quantization of fields in dielectrics to spontaneous 
emission is that the oscillators that comprise the dielec- 
tric host are assumed to be unaffected by the presence 
of the embedded atom. The dielectric medium, as well 
as the field, is treated as a local condition at the site 
of a resonant atom such that the atom interacts with 
an all-pervasive, nonlocal, quantized effective field, the 
vacuum polariton modes, rather than the local vacuum 
radiation field modes and the oscillators. Because the 
near-dipole-dipole interaction is the fundamental inter- 
action underlying the Lorentz local field, a many-body 
approach that explicitly deals with the vacuum radiation 
field modes and the interactions of the atom with the 
nearby polarizable particles of the host is clearly needed 
to accurately evaluate the effects of local fields on spon- 
taneous emission in dielectric media. 

In this letter, we develop Langevin-Bloch operator 
equations of motion for a dense collection of two-level 
atoms embedded in a dielectric host medium. The 
Heisenberg picture is used since this provides the most di- 
rect correspondence between the operator equations and 
the Bloch equations. We begin the development from a 
fully microscopic many-body viewpoint in which the ma- 
terial is treated as a disordered mixture of two different 
species of two-level systems and derive the Heisenburg 
equations of motion for the material and field mode op- 
erators. Adiabatically eliminating the variables associ- 
ated with the quantized field modes results in coupled 
equations of motion for the material variables. We take 
the harmonic oscillator limit for one species by assum- 
ing that its resonance frequency is sufficiently detuned 
from the primary species that the atoms remain in the 
ground state. Adiabatically eliminating the harmonic os- 
cillators results in a Langevin-Bloch formulation for two- 
level systems embedded in a dielectric host that exhibits 
local-field renormalization of the fluctuations, the near- 
dipole-dipole (NDD) interaction of the two-level atoms, 
the radiation field, the dephasing rate, and the popula- 
tion decay rate. We obtain 
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for the renormalized spontaneous emission rate. In our 
many-body derivation, the dielectric local-field enhance- 
ment factor £ = (n 2 + 2)/3 arises from the interaction 
of the embedded atom with the nearby polarizable parti- 
cles of the host via the electromagnetic field. The linear 
index of refraction is complex and frequency dependent 
and properly accounts for the frequency dispersion and 
absorption of the dielectric. 
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We consider a disordered mixture of two species, a and 
6, of two-level systems. The two-level systems are cou- 
pled only via the electromagnetic field. We allow for the 
possibility of an externally applied probe or driving field 
that is taken to be in a coherent state. The constituents 
of the total Hamiltonian are the Hamiltonians for the 
free atoms of species a and b, the free quantized radi- 
ation field, and the interaction of the two-level systems 
with the quantized electromagnetic field. We have, in the 
electric-dipole and rotating-wave approximations, 

j n l,a 
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where a\ and a; are the creation and destruction op- 
erators for the field modes and uji is the frequency of 
the field in the mode I. The vacuum dispersion rela- 
tion is used throughout, e.g. ki = kiuii/c, where k\ 
is a unit vector in the direction of fc;. For atoms of 
species a, <j\ is the population inversion operator and crj_ 
are the raising and lowering operators for the j atom, 
gj = (2TTU>i/hV) 1 ^ 2 ^ a f)j-ej: i a is the coupling between the 
atom at position fj and the quantized radiation field, pj 
is a unit vector in the direction of the dipole moment at 
fj, io a is the transition frequency, fi a is the dipole mo- 
ment, il a — Ha£/h is the Rabi rate, and £ is the field 
envelope associated with the coherent field with carrier 
frequency u> p . For species b, Q, (±, hf, r„, p„, u> b , 
and fib — [ib£/h respectively perform the same functions. 
Finally, V is the quantization volume, eg g is the polar- 
ization vector, and a denotes the state of polarization. 

The Heisenberg equations of motion for the material 
and field mode operators 



are derived from ih(dO/dt) — [0,H]. From this point, 
we adopt normal ordering in which a\ appears to the left 
of the atomic operators and a; appears to the right. 

The Heisenberg equations of motion for the material 
variables are obtained by elimination of the variables as- 
sociated with the quantized field modes [fill. Then, 
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Using standard QED methods in the Markovian approx- 
imation |12-n3], this immediately reduces to 
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in a frame of reference rotating at tu p , such that A a = 
lu p — uj a . In the context of the analysis of Ref. Jl2| , 
with normal ordering, / + is a Langevin force operator 
arising from fluctuations of the vacuum field, the de- 
phasing rate 7 a /2 is half of the population decay rate 
a{^t a e~ i{ - Upt ~ kp ' ? ^ la = ^p\fJ-a\ 2 /(Sc 3 h) that is asssociated with the self- 
field under the condition i = j, and e a — 4:irN a \/j, a \ 2 / (3h) 
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is the strength of the near-dipole-dipole (NDD) interac- 
tion due to the reaction field of all other atoms, i ^ j, of 
species a, where N a is the relevant number density and 
<r_ is a local spatial average of the operator JlJ,|l^] . The 
last term in Eq. (^) is the contribution of the reaction 
field from all the atoms of species b to the j th atom of 
species a. A similar analysis yields 
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for the equation of motion of the inversion operator. 
In the harmonic oscillator limit for species b, we have 
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47r7Vfc|/Lt{,| 2 /(3?i). Equations fl3|)-(|5D are coupled opera- 
tor equations for a material composed of two-level sys- 
tems and harmonic oscillators. None of the relevant pa- 
rameters for the two-level systems, i.e. frequency, field 
strength, fluctuations, dephasing rate, population decay 
rate and NDD, are renormalized by the prescence of the 
host medium as long as we retain separate equations of 
motion for the harmonic oscillators. The next step is to 
adiabatically eliminate the equations of motion for the 
harmonic oscillators by substituting the formal integral 
of Eq. (g) into Eqs. © and (|). Thus 



of the original field components, the coherent field, vac- 
uum fluctuations, the self-field, and the reaction field, 
and will lead to the renormalization of each. The self- 
field contribution, i — j, in which the atom couples to 
itself via the linear particles, can be evaluated using the 
transverse delta function || . The reaction field contribu- 
tion, the near dipole-dipole interaction of all the i atoms 
with atom m, is of the same form as the NDD interac- 
tion of the i atoms with atom j that was studied in Ref. 
]l4j and used in obtaining Eq. (g). We refer to this type 
of interaction as a many-atom Milonni-Knight problem 
[O.E5|. The sum in the large square brackets becomes 



In the last term of Eq. (|J), note that Jl3| 
dt"e^'-^f(t") -> —/(*'), 
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since the exponential is strongly peaked near t' = t" . 
Because r m ^ fj, the remaining part of the last term in 
Eq. (||) is again the many-atom Milonni-Knight problem. 
Then the portion exterior to the large square brackets can 
be written as 
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where N b is the number density of oscillators, species b 
. nl + 2 



(10) 



is the complex Lorentz dielectric local-field enhancement 
factor, and n b is the complex index of refraction of the 
host material. Using (^) and (^) in Eq. (^), we obtain 
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A similar calculation for the equation of motion of the 



where £+ (0) = and a = i(A b + e b + ry 6 /2). 

The last term is the part of the reaction field that inversion operator, Eq. (^), yields 
is due to the prescence of the harmonic oscillators and 
their subsequent adiabatic elimination. The large square 
bracket contains terms that are largely equivalent to all 



dt 



3 



-l la {ai + l) + H.c. 



(12) 



Equations @ and (|T|) can be considered as operator 
Langevin-Bloch equations of motion for a dense collec- 
tion of two- level atoms embedded in a dielectric medium. 
The effect of the adiabatically eliminated damped linear 
oscillators is contained in the complex Lorentz dielec- 
tric enhancement factor I that renormalizes the coherent 
field, the Langevin force operator, the NDD interaction, 
the dephasing rate and the population decay rate. We 
have neglected Cauchy principle parts throughout and 
there will be local-field enhancement effects from these, 
as well, e.g. renormalization of the Lamb shift by Re(£). 

The results of our fully microscopic many-body QED 
treatment agree with semiclassical results for local-field 
enhancement of the coherent field Jlj] and the NDD inter- 
action , if expectation values are taken in the limit of 
classical factorization. The purely quantum mechanical 
aspects are the enhancement of the Langevin force oper- 
ator and the corresponding damping rates. At resonance, 
the population decay rate that appears in the Langevin- 
Bloch equation of motion for the inversion operator is the 
renormalized spontaneous emission rate 
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The spontaneous emission rate is renormalized by Re(£). 
In addition, there is a level shift, or frequency renormal- 
ization, in the amount of Im(£)"f a /2 due to the product 
of the imaginary part of I with the dephasing rate. 

There have been a number of experimental mea- 
surements of the spontaneous emission rate of embed- 
ded atoms, or atom-like particles, in a dielectric |[L7) . 
However, these experiments typically involve non-trivial 
boundary conditions, such as organic ligand cages or 
nanospheres, that can profoundly affect the observed 
spontaneous emission rate. To date, we know of no ex- 
perimental measurements of the index dependence of the 
spontaneous emission rate in a bulk dielectric. The com- 
plete theory presented here makes it possible to bring 
the entire arsenal of laser spectroscopy to bear on the 
measurement of local-field effects due to a dielectric back- 
ground. For example, because even small frequency shifts 
can be resolved spectroscopically, it might be possible to 
verify our results by measuring the level shift Im(£)"/ a /2 
or the renormalization of the Lamb shift by Re(£) as a 
function of the density of a buffer gas. 

Operator equations of motion, including the local-field 
renormalization of the spontaneous emission rate, for a 
single atom, or a tenuous collection of atoms, embed- 
ded in a dielectric are contained in these, more general, 
results in the limit e a — > 0. The NDD interaction of 
densely embedded atoms in a dielectric is included pri- 
marily because it provides a very useful theoretical back- 
drop for evaluating the many-body effects. The same 



type of many-body summations that occur in the local- 
field renormalization of the spontaneous emission rate 
have been studied previously in the evaluation of local- 
field effects for dense atoms in vacuo jl4|. Further, it was 
our semiclassical derivation of the 'anomalous' renormal- 
ization of NDD interaction in a dielectric host jl(J that 
indicated a need to examine the problem in its entirety. 
Because the NDD interaction and the spontaneous emis- 
sion rate have the same dependence on the dipole mo- 
ment and arise in the same way from the elimination 
of the field operator, one would expect them to have the 
same renormalization in a dielectric. We have shown that 
this is the case. 

In conclusion, we obtained the renormalization of the 
spontaneous emission rate of an atom embedded in a di- 
electric host. This result was obtained from a fully mi- 
croscopic, many-body derivation of Langevin-Bloch op- 
erator equations for two-level atoms embedded in an ab- 
sorptive and dispersive, linear dielectric host. The di- 
electric local-field enhancement of the coherent field, the 
Langevin force operator, the NDD interaction, and the 
damping rates all stem from the same reaction field that 
arises from the nearby harmonic oscillators, necessitat- 
ing the full many-body derivation. We found that the 
dielectric enhancement of the spontaneous emission rate 
is much smaller than indicated by previous studies. This 
is an enabling result, paving the way for application of 
high-index materials to enhance quantum optical effects. 
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